Persistent random walk on a site-disordered one-dimensional lattice: Photon 

subdiffusion 
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We study the persistent random walk of photons on a one-dimensional lattice of random transmit- 
tances. Transmittances at different sites are assumed independent, distributed according to a given 
probability density f(t). Depending on the behavior of f(t) near t = 0, diffusive and subdiffusivc 
transports are predicted by the disorder expansion of the mean square-displacement and the effec- 
tive medium approximation. Monte Carlo simulations confirm the anomalous diffusion of photons. 
To observe photon subdiffusion experimentally, we suggest a dielectric film stack for realization of 
a distribution f(t). 
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I. INTRODUCTION 

Random walks figure prominently in a multitude of 
different physical problems. This is exemplified by such 
diverse fields as the polymer physics |I| , crystallographic 
statistics |2j, |3( , transport in disordered media |3|, U, |5j , 
bacterial motion, and other types of biological migration 

@ 

Random walks with correlated displacements, model a 
host of phenomena. For example, the vacancy mechanism 
of atom diffusion in solids incorporates a correlation ef- 
fect, since an atom has a larger probability to move back- 
ward to the hole it just vacated rather than onward 0. 
Correlations also arise in diffusion of guest molecules in 
zeolite channels Q, electron hopping in Coulomb glass 
0, motion of excitons at low temperatures in mixed 
naphthalene crystals ^(j; e t c - Among the correlated 
walks, the persistent random walk is possibly the sim- 
plest one to incorporate a form of momentum in addi- 
tion to random motion |3|. . In its basic realization on a 
one-dimensional lattice, a persistent random walker pos- 
sesses constant probabilities for either taking a step in 
the same direction as the immediately preceding one or 
for reversing its motion. Firstly introduced by Furth as 
a model for diffusion in a number of biological problems 
[Tl| , and shortly after by Taylor in the analysis of turbu- 
lent diffusion [lj], the persistent random walk model is 
now generalized to study, e.g., polymers 0], chemotaxis 
fl4|. cell movement ^5j, an d general transport mecha- 
nisms muni- 

Recently, the persistent random walk model is utilized 
in the description of diffusive light transport in foams 
l2ll 122^ which is well established by experi- 



ments [23|. A relatively dr 
rated by thin liquid films 



foam consists of cells sepa- 
. Cells in a foam are much 
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larger than the wavelength of light, thus one can employ 
ray optics and follow a light beam or photon as it is trans- 
mitted through the liquid films with a probability t called 
the intensity transmittance. This naturally leads to a 
persistent random walk of the photons. Special attention 
is paid to the light transport in the ordered honeycomb 
and Kelvin structures, which have been used for an ana- 
lytic access to the physical properties of disordered foams 
as exemplified by work on their rheological behavior [25| . 
Symmetries of hexagonal and tetrakaidecahcdral cells al- 
low identifying specific one- and two-dimensional random 
walks, which are absent in disordered foams 0, l2lj . 
Moreover, in a first model, it is assumed that intensity 
transmittance does not depend on the incidence angle 
of photons and film thickness. Nevertheless, analytical 
treatment of these peculiar and simple random walks 
facilitates interpretation of Monte Carlo simulations, in 
which topological disorder of foams and exact thin-film 
intensity transmittance are taken into account |19l |21| . 

In the ordered honeycomb (Kelvin) structure, the one- 
dimensional persistent walk arises when the photons 
move perpendicular to a cell edge (face) . Thin-film trans- 
mittance depends on the film thickness. Films are not 
expected to have the same thickness. These observa- 
tions motivates us to consider persistent random walk 
on a one-dimensional lattice of random transmittances. 
We assume that transmittances at different sites are in- 
dependent random variables, distributed according to a 
given probability density f(t). 

Our first approach to the problem is based upon a dis- 
order expansion of the mean square-displacement due 
to Kundu, Parris, and Phillips [2(|. Assuming that 

< 1/t >= L f(t)/t dt is finite, we validate the classi- 
cal persistent random walk with an effective transmit- 
tance t e ff, where l/t e // =< 1/i >■ Inspired by this 
result, we generalize the effective medium approxima- 
tion (EMA) formulated by Sahimi, Hughes, Scriven, and 
Davis [53, [28] to investigate the transport on a line with 
infinite < 1/t >. We show that if f(t) -> /(0) as t -> 0, 
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the mean square-displacement after n steps is propor- 
tional to n/ln(n). If f(t) ~ f a t~ a (0 < a < 1) as 
t — > 0, we find that the mean square-displacement is pro- 
portional to n ( 2 - 2a )/( 2 - Q ). Our Monte Carlo simulations 
confirm the anomalous diffusion of photons. Quite in- 
teresting, we find that anomalous diffusion of persistent 
walkers and hopping particles on a site-disordered lattice 
l28l |2jj are similar. Finally, for the experimental 
observation of the photon subdiffusion, we suggest a di- 
electric film stack to realize small transmittances. 

Our paper is organized as follows. In Sec. [H] we in- 
troduce the model. The perturbative approach and the 
effective medium approximation are discussed in Secs. lIIll 
and llVI respectively. The numerical treatment and its re- 
sults are reported in Sec. [V] We close with a discussion 
of our results and conclusions in Sec. IVII 



II. MODEL 

We consider a one-dimensional lattice random walk in 
which steps are permitted to the nearest neigbor sites 
only. We normalize the length and duration of a step 
to I. At each site j, a walker either takes a step in the 
same direction as the immediately preceding one with a 
probability tj, or reverses its motion with a probability 
r j = 1 — tj. Here we assume symmetric transmittances, 
i.e. tj^j+i = tj^j-i = tj, as the transmittance of a 
thin-film is the same whether the light ray is going to 
the right (+) or to the left (— ) direction. 

We assume that i) transmittance at each site is a ran- 
dom variable, ii) transmittances at two different sites 
are independent, iii) transmittances at all sites are dis- 
tributed according to a given normalized probability den- 
sity f(t). Apparently J Q f(t)dt = 1. For any function 
h(t), we define < h(t) >= h{t)f{t)dt. 

We denote by P + (n, j) (P~(n, ji)) the probability that 
the walker after its nth step arrives at site j with positive 
(negative) momentum. A set of two master equations can 
be established to couple the probabilities at step n + 1 to 
the probabilities at step n: 

P+(n + l,j) = tj^P+(n,j - I) + rj -iP-(n,j - I), 
p-(n + l,j) = rj +1 P + {n,j + l)+t j+1 p-(n,j + l). 

(1) 

For the description of the photon distribution on the 
line, we do not need to specify the internal state (±) 
explicitly. That means we are mainly interested in the 
probability that the photon arrives at position j at step 
n, 

P(nJ)=P+(n,j)+p-(nJ), (2) 
from which we extract the first and second moments after 



n steps as the characteristic features of a random walk: 

«J»n = <E-7' P ("> ■?')>> 
3 

((3 2 ))n = C£j 2p (nJ))- (3) 
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Here the first bracket represents an ensemble average over 
all random transmittances, and the second bracket signi- 
fies an average with respect to the distribution P(n,j). 

One obtains the classical persistent random walk as- 
suming a constant transmittance t at each site. Trans- 
lational invariance of the medium is then invoked to de- 
duce the exact solution of P(n,j) in the framework of 
characteristic functions (the spatial Fourier transforms of 
probability distributions) 0, y| . Furthermore, the mean 
square-displacement of photons after n — > oo steps can 
be obtained as 

(f)n = -^-n. (4) 

Considering a lattice of random transmittances with 
an almost narrow distribution f(t), one may intuitively 
expect normal diffusion of the photons, and the valid- 
ity of Eq. Q with an effective transmittance < t >. 
However, a closer inspection reveals that even a few sites 
with small transmittances (large reflectances) may dras- 
tically hinder the photon diffusion: In the extreme limit 
where two transmittances are zero, photons are confined 
between two sites. This peculiar aspect of diffusion on 
a one-dimensional lattice, rules out the above guess. In 
the following section, we present a sound perturbative 
approach to the problem. 

III. SYSTEMATIC DISORDER EXPANSION 
OF MEAN SQUARE-DISPLACEMENT OF THE 
PHOTONS 

Many of the approaches to the transport in disordered 
media have the disadvantage of being restricted to one- 
dimensional problems. Here we adopt the method of 
Kundu et al. [2(| , which is applicable to two- and three- 
dimensional media. We obtain a series of approximate 
solutions for photon transport on a disordered line by 
transforming the master equation Q to an equivalent 
but more appropriate integral equation for the charac- 
teristic function 

oo 

P ± (n,0)= P ± {n,jY je . (5) 

j—-oo 

This crucial step can be achieved utilizing the generalized 
generating functions 

P ± (n,6)=J2tjP ± (^3)e lje , (6) 

3 

as will be shown in the following. 
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First we simplify the set of coupled linear difference Thus P ± (z,0) can be easily expressed in terms of 

equations QJ using the method of the z-transform [U|3(j P ± (z, 0). From Eqs. JSJ and © 
explained in Appendix 

P+(z,j) P+(n = 0,j 



P-{z,j) P-(n = 0,j) 



= tj^P+izJ-l) 
+r j - 1 p-(z,j-l), 

= r j+1 P+(z,j + l) 

+t j+1 p-(z,j + l). (7) 



We assume the initial conditions P + (n = 0, j) = P~ (n — 
0>j) = <5j,o/2. Now we multiply both sides of J7J) with 
exp(ij8) and sum over j, which leads to 



M 2 (6>) 



P+(z,0) 

p-(z,e) 



•Mi(«), 



(8) 



where 



M 3 (z,0) 
M 2 (0) 



- — e 
l 



P±(z, 0) = £ -i- f P±(z, 0)e-W- e )d0, (10) 



hence we can immediately transform Eq. (|SJ) to an equiv- 
(9) alent integral equation for P ± (z, 0). It reads 



« I > I- M 3 -.( 2 ,.,M 2W ) (|!<^)) . M 3 -(M)M lW - ? ££ e — »> (|!< :; *>) * 



(11) 



where the random variable Aj is defined as 



A, 



1 

< - >, 



(12) 



and I is the identity matrix. Note that two assumptions, 
that none of the transmittances is zero, and < 1/t > is 
finite, are required for the validity of Eqs. ijTU|) - ((T2"|l . 

Successive approximations to the solution of integral 
equation (|llfl can be generated by the iteration method. 
Then ensemble averaging over the random transmit- 
tances is straightforward. We use the identity < Aj >= 
and the assumption 



(13) 



where matrix M^z, 6) is given in Appendix iBl 
P o (z,0) = 



1 z{< \ > e w - 2cos0)+ < \ > 



2 < i > z 2 (2- < \ >) - 2zcos6>+ < \ >' 

(15) 

and * denotes the complex conjugation. 

The task is now calculating < P ± (z,9) > from Eq. 
©, and using the identity 

EW..I) = l^^\e=° = 1.2,...). (16) 

3 

to obtain mean square-displacements of the photons. To 
this end, 



£«j 2 »n* T ' 



z - 4iz( < 



dP+(z,0) dP-(z,0) _ \, 

— Be W~ > )\ e =° 



to obtain the leading terms: 



<(t + (z> e A >= (l + A»M,(M)) (^ 6) X (14) 



+ 



< i >-i(l- zf 



V2A 2 



2(< \ > -l) 2 - 5 (1-z) 1 - 5 ' 



(17) 



4 



((f)) 



1. We use Tauberian theorem l|A2(l to conclude 
1 V2A 2 



<7>-l 



V^(< 7 > -I) 2 



.5 



(18) 



Mean square-displacement of the photons after n — > oo 
steps is indeed proportional to n. We validate the clas- 
sical persistent random walk result with an effective 
transmittance t e f f , where 



11 

teff ^ 



(19) 



To facilitate solution of Eq. (JTJ) we introduce a refer- 
ence lattice or average medium, with all intensity trans- 
mittances (reflectances) equal to t e {z) (r e (z)), and prob- 
abilities P e ± (z, j), so that 



P+(z,j) P+(n = 0J) 



P-(z,j) P-(n = 0,j) 



t e (z)P+(z,j-l) 
+r e (z)p-(z,j-l), 
r e (z)P+(z,j + l) 
+t e (z)p-(z,j + 1)(20) 



IV. EFFECTIVE MEDIUM APPROXIMATION 

Stochastic transport in random media is often subd- 
iffusive or superdiffusive [H 0> IH EH ■ At this stage, we 
pay attention to persistent random walk on a line with 
infinite < 1/t >, where the normal diffusion is not guar- 
anteed. Our approach to the problem is based upon a 
variant of effective medium approximation (EMA) devel- 
oped in H^lp. 



EMA determines t e (z) and r e (z) in a self-consistent man- 
ner, in which the role of distribution f(t) is manifest. 
This is done by taking a cluster of random transmittances 
from the original distribution, and embedding it into the 
effective medium. We then require that average of site 
occupation probabilities of the decorated medium dupli- 
cate P^(z, j) of the effective medium. We will sketch the 
method in the following. 

Subtracting Eqs. Q and i|2U[l . we obtain 



1 



Q+(z,j) 

t 



'3-1 





j-1 





T~(z) 

-T-(z)] 



Q + (z,j-l) 
Q-(z,j-l) 

P+{z,j-l 
P~(z,j-1 



-T+(*) 



} + (z,j + l 
)-(z,j + l 





3+1 



-T+(z)] 



p+(zj + r 

P-(z,j + l] 



(21) 



where 



) + (z,j) 

T-(z) 
T+(z) 



[P + (zJ] 
\P-(z,j) 

t e {z) r e {z) 




r e (z) t e (z) 



Pe(zJ) 



(22) 



Equation 121|) suggests to define an associated Green 

Ull (jrl2 
G*21 G22, 



function G(z,j) 



by the equation 



-G(z,j) - T-(z)G(z,j - 1) - T+(z)G(z, j + 1) 

z 



whose solution is the inverse Fourier transform of 

~2 



For the present, we consider only the simplest approxi- 
mation, and embed one random transmittance at site I 
of the effective medium. Then solution of Eq. (|21|l is 



Q (z,j] 



where 

S(z,0) = 



G(z,0)S(z,8)e- wu - l) 



P+(z,l)\ dd 



P-(z,l)l 2ir 



[U-UWe" [n-r e (z)]e id 
n - r e (z)]e-* s [i, - t e (z)]e-* e 



KZ> ( 25 ) 



(26) 



^'. oI 'or 
(23) 



Self-consistency equation is < P ± (z 7 l) >= P^(z,l) 



/■2tt jrt 

<[i-y o G(z,e)s(z,e)^-]- 1 >=i. 



(27) 



G(z,6) 



1 - 2zt e (z) cos 6 + z 2 [tl (z)-rl{z)] 
' x -t e {z)e' ie r e (z)e l9 
r e (z)e-* e \-t e {z)e ie 



(24) 



The above matrix equation leads to two independent self- 
consistency conditions. Choosing t e (z) + r e (z) = 1, one 
of the conditions can be satisfied. It is a signature of 
the success of effective medium to represent the original 
random medium, where tj + r j = 1 holds at any site 
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j. The second self-consistency condition then determines 

t e (z): 



1 



f(t)dt 



te(z) Jo [t e ( Z )-t]-. 



y/l-[2t e (z)-l]2z- 



+ t 



(28) 



in which the role of distribution f(t) is manifest. 

Z-transform of mean square-displacement of the pho- 
tons in the effective medium can be obtained as 



oo 

E 

n=0 



«i 2 »»* n = 



z 1 + z[2t e (z) - 1] 
(l-z) 2 l~z[2t e (z)-l]' 



(29) 



We are interested in the long time behavior, thus Taube- 
rian theorems suggest to analyze Eqs. I|28|l and l|29|l in 
the limit z — ► 1. 

First we assume that < l/t > is finite and t e (z) has 
no singularity in the limit z — > 1. Then Eq. I|28|) yields 



t e (*) 



1 /(t)dt 



(30) 



in accordance with our second assumption. We deduce 
from Eqs. J2HJ) and (|X2|) that 



«J 2 »n = 



<7>-l 



n. 



(31) 



therefore the system evolves diffusively. This is in com- 
plete agreement with the predictions of Sec. IIIII 

When < l/t > is infinite, the behavior of l/t e (z) is 
determined by the behavior of f(t) at small values of t. 
Let us assume that f(t) has a finite derivative at t = 0. 
We can decompose the integral in Eq. (|28Jl into a sum: 



f(0)dt 



te(z) Jo [U Z )-t]-, +t 

1 ^/l-[2t c (z)-l] 2 z 2 

[f(t) - f(0)]dt 



o \Uz)-t}—, — £p£ +t 

1 CV 1 1 y/l-[2t e (z)-l]*z2 



f(t)dt 



y/l-[2U(z)-lYz- 



+ t 



-, (32) 



where e is a small number. By assumption f(t) — /(0) = 
tf'(0) near t — 0, and the factor t cancels the potential 
singularity in the second term as z — * 1. Indeed the only 
singular behavior can come from the first term. Equation 
()32|1 then leads to 



t e (z) 



-/(0)ln[ 



t e (z)VT~^ 



sfl - [2t e ( Z ) - 1]2 Z 2 



(33) 



Methods for the asymptotic solution of transcendental 
equations |32| are invoked to obtain 



/(0) ln(l-z) 



(34) 



We deduce from Eqs. (HHJ) and (|A3|> that 

2 n 



((3 2 ))n = 



/(0) ln(rc)' 



(35) 



thus the transport is subdiffusive. 

We now investigate the cases in which /'(0) is infinite. 
If f(t) ~ f a t~ a (0 < a < 1) as t — ► 0, self-consistency 
equation l)28|) yields 



[*e(*)] 



a-1 



7r/ a r l-[2i e (z)-l] 2 z 2 lf 



Sm(7TQ!) 



1-Z 2 



(36) 



where we have used Jg 00 x a /(l + a^)da; = 7r/ sin(7ra). In 
this case we find 

^)-[2-/ 2 ^]^(l-,)^, (37) 



and 



(0' 2 ))n 



- n 2 -° 



p/ 4— jo \ 



(38) 



Apparently, < (2 — 2a)/(2 — a) < 1 and the transport 
is subdiffusive. 



V. NUMERICAL SIMULATIONS 

The predictions of EMA can be inspected by numer- 
ical simulations. The computer program produces 400 
media, whose transmittances are distributed according 
to a given f(t). For each medium, it takes 10 3 photons 
at the initial position j = and generates the trajec- 
tory of each photon following a standard Monte Carlo 
procedure. 

For the binary distribution f(t) — p\8[t — t\) + (1 — 
pi)5{t-l + ti), where pi £ [0.1, 0.2, 0.9] and t x £ 
[0.1,0.2,0.3,0.4], statistics of the photon cloud is eval- 
uated at times n £ [7000, 7100, 10000]. The mean- 
square displacement ((j 2 )) n is computed for each snap- 
shot at time n, and then fitted to Dn + O by the 
method of linear regression. An offset O takes into 
account the initial ballistic regime. Fig. ^ shows the 
excellent agreement between numerical simulations and 
D(pi,ti) = ti(l - h)/( Pl - 2pih + <?), predicted by Eq. 
(|31|l . The maximum differences ±0.04 are comparable to 
the errorbars ±0.01 which linear regression anticipates. 

For the uniform distribution f(t) = 1 (0 ^ 
t 1), statistics is evaluated at times n € 

[10000, 25000,. ..,400000]. In Fig. ((j 2 »„/n, and 
2/ln(n) are plotted. This figure confirms the EMA pre- 
diction 1351) . 

For f(t) = (1 - a)t- a (0 < t < 1), where a £ 
[0.1, 0.2, 0.9], statistics is evaluated at times n £ 
[70000, 71000,. ..,100000]. ln[((j 2 ))„] is computed for 
each snapshot at time n, and then fitted to 7ln(n) + O, 
where O is an offset. In Fig. [3 we plot 7 and (2 — 2a)/(2— 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Pi 

FIG. 1: The diffusion constant as a function of pi and ti 
which parametrize f(t) = p\S(t — t\) + (1 — pi)S(t — 1 + ti). 
Theoretical and Monte Carlo simulation results are denoted, 
respectively, by lines and points. 




0.15 



10000 100000 190000 280000 370000 
n 

FIG. 2: {{f)) n /n as function of n, for fit) = 1 (0 < t < 1). 
Theoretical and simulation results are denoted, respectively, 
by line and points. 



a) as a function of a, to inspect the EM A prediction (|38ll . 
The differences are greater than the errorbars which re- 
gression anticipates, but the overall agreement is excel- 
lent. 



VI. DISCUSSIONS AND CONCLUSIONS 

In the present paper, we address the persistent ran- 
dom walk on a one-dimensional lattice of random trans- 
mittances. The photon transport is diffusive, provided 
that < 1/t > is finite (class I). The transmittance of the 
effective medium is given by Eq. I|I9|) . It expresses the 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

a 



FIG. 3: 7 = hi[((j 2 ))n]/ ln[n] as n — » oo, versus a which 
parametrizes f(t) = (1 — a)t~ a . Theoretical and simulation 
results are denoted, respectively, by line and points. 



fact that a few sites with small transmittances (large re- 
flectances) may hinder the photon diffusion. As percola- 
tion properties 28] , this feature is induced by the dimen- 
sionality of the lattice. Furthermore, Eq. (|I8|I shows that 
fluctuations in transmittances give rise to a ^/n correc- 
tion to the mean square-displacement. Nevertheless, this 
correction is negligible at times n ^> A 4 /(< 1/t > —I) 3 . 

A photon steps back by each reflection. Intuitively, one 
expects the abundance of large reflectances to drastically 
decrease excursion of the photons. Our EMA predicts 
a subdiffusive transport when < 1/t > is infinite. The 
self-consistency equation (I32[l divides such distributions 
to two classes II and III, which are summarized in Table 
[I] Remarkably, for class III the exponent of mean square- 
displacement is distribution dependent (nonuniversality) . 

It would be instructive to compare transport of per- 
sistent walkers and excitation dynamics in random one- 
dimensional systems |53, IS . Hopping conduction is 
described by the master equation 



dP(T,j) 



= Wj j^[P(T,j-l)-P(T,j)} 



+ Wj . j + 1 [P(T,j + l)-P(T,j)}, (39) 

where ^ w j.j+i — w j+i.j ^ oo is the random transition 
rate between sites j and j+1, and r is the continous time. 
Correlated steps are the essence of Q , while independent 
steps are the base of Despite this celebrated differ- 
ence, Table 13 reveals remarkable similarities. 

As we mentioned in SecQJ the persistent walk on a one- 
dimensional lattice of random transmittances arises when 
the photons move perpendicular to an edge (a face) of the 
honeycomb (Kelvin) structure. At normal incidence on 
a film of random thickness dyy and refractive index nw. 



t = 1 



2r^(I - cos p w ) 



w 



2r^ cos (3 W 



(40) 
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TABLE I: EM A prediction of ultimate growth with time of the mean square-displacement. 



Photons 


Excitons a 


/(*) (U 2 ))n 

I < j oo n 

II < j >= oo, /(£) -> /(0) as t -> ra/ln(n) 

III /(f) ~ f- Q as t n (2-3a)/(2-a) 


/H ((3 2 ))r 

< i >= oo, /(w) -> /(0) as iu -»• r/lnr 

/H ~ lo"" as W -» r (2-2a)/(2-a) 



a From 27]. r and u> denote the continous time and the transition rate between the nearest neighbours, respectively. 



where rw = {nw — + l)n (3w = ^irdwnw I A, and 

A is the light wavelenght |L9t |21|. For foams nw ~ 1.34 
(water) and rw ^ ±1, hence the above transmittancc 
never approaches zero, < 1/t > exists, and the photon 
transport is diffusive. 

In a realization of a distribution /(f) ~ f a t~ a as f — * 0, 
small transmittances prevail. For an experimental obser- 
vation of the photon subdiffusion, we propose dielectric 
mirrors to obtain desired small transmittances. Consider 
the multilayer cofiguration 



HL HL- 



N 

HL W LA- 



iY 

LH LH. 
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APPENDIX A: z-TRANSFORM 

The z-transform F(z) of a function F(n) of a discrete 
variable n — 0,1,2, ... is defined by 



where H (L) is a quarter-wave layer with high (low) index 
n>H and N is the number of HL (LH) pairs. W is 

a layer of thickness dw and index nw- Following the 
theory of multilayer films |33j |. we obtain 



t = 1- 



(rV-l) 2 (l-cos/3 w ) 



r' w {l — cos/%) + 2r' w (3 + cosflw) + 1 - cos [3w 

(41) 

where r'w — n w( n H/ n L) AN ■ Apparently, as N increases 
the transmittance rapidly approaches zero. For exam- 
ple, suppose n H = 2.40 (Ti0 2 ), n L = 1.38 (MgF 2 ), 
nw — 1-46 (Si02), A ~ 550nm and fiw — t (quarter- 
wave layer). For 1, 2 and 5 pairs, the transmittances are 
0.185604, 0.022174 and 0.000029, respectively. In this re- 
alization of small transmittances, two points are consid- 
ered. First, transmittance of the multilayer is the same 
whether light rays go to the right or the left direction. 
Second, parameters of the multilayer, especially N and 
dw, can be varied to fine tune the transmittance. More- 
over, absorption of the multilayer is negligible. Note that 
thin metalic slabs are highly reflecting, but absorptive. 

Our studies can be extended to higher dimensional lat- 
tices, and special two-dimensional photon paths in the 
honeycomb and Kelvin structures 0, . Another path 
to pursue is the creation of artificial one-dimensional 
structures to observe anomalous diffusion of photons. 
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F{z) = Y j F{n)z n . 



(Al) 



n=0 



One then derives the z-transform of F(n + 1) simply as 
F(z)/z — F(n = 0)/z. Note the similarities of this rule 
with the Laplace transform of the time derivative of a 
continuous function 0, |3(j • 

Under specified conditions the singular behavior of 
F(z) can be used to determine the asymptotic behav- 
ior of F(n) for large n (Tauberian theorems) 0|- For 
example: 



F(z) 



F(z) 



r(i-q) 
i 



ln( T L-)(l-z)« 
where r(a) = J °° e~ t t a - 1 dt. 



F{n) 
F(n) 



— ( A2 ) 



r(a) ln(n) 



(A3) 



APPENDIX B 



The matrix ~M.n(z,6) introduced in Eq. (|14|> . can be 

conveniently written as 



M 4 (z,0) = 



a(z,0) + b(z,6) a(z,6) -b(z,0) 
a(z,6)-b*(z,6) a(z,0) + b*(z,d) 



(Bl) 
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where 

a(z,d) = 

b(z,0) = 




-z 2 (< i > -2) 2 
z 2 < \ > -2?zsni0- < \ > 
z 2 (2- < \ >) -2zcos0+ <\>' 
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